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1. Introduction 

The K —> 3n and the v\ — > 3ti decays are important sources of information for the studies of 
isospin breaking effects. From the branching ratio of the v\ — > 3% decay, which is forbidden in the 
isospin limit, we can obtain values of the isospin breaking parameters like R = ™ s d -Z ■ The final 
state interactions of two neutral pions coming from the K — > 3n or the v\ — > 3ti decays resulting 
in the so-called cusp effect offer quite a simple experimental determination of the 5-wave 7171 
scattering lengths. This has been already successfully done for the K + — > 7i + n°7l° decay and 
nowadays there appear a large effort in a careful measurement of the cusp effect also in the case of 
K l ^3k° and 77 — > 3tt° decays. 

These are just few of the reasons why these decays attract so much attention and are still under 
intensive studies, both from the experimental and from the theoretical points of view. For the sake 
of limited space here, we refer the reader for a more comprehensive list of these works to our 
previous proceedings Nevertheless, for the reader's convenience we list here the experimental 
collaborations undertaking these measurements |2|] and theoretical studies dealing with T/ — ^ 3?r 
decay [Q]-[|5|] as we concentrate here mainly on results of this decay. 

The processes in question are (generically denoted as P — > 7C7r7r in the following section) 

K + -» 7r+7r 7r°, 7r+7r + 7r"; K L -> 7r 7r°7r°, 7r°7z; + 7zr; tj -» 71°k°k°, 7T°7z; + 7r; K s ^n°n + n-. 

(1.1) 

The standardly used model-independent parametrization of the energetic dependence of the 
Dalitz plot for these decays (Dalitz parametrization of PDG) does not account for the analytic 
structure arising from final state interactions in the presence of isospin breaking. 

2. Parametrization and reconstruction procedure 

We propose a unified model-independent parametrization of the corresponding amplitudes 
valid at two loops relying on very general principles, unitarity, analyticity and crossing symmetry, 
combined with chiral counting. Our treatment is fully relativistic and the 7T7T scattering parame- 
ters appearing there can be chosen so that instead of their chiral expansion, we use their physical 
interpretation (values) up to (and including) two-loop order. 

The parametrization takes the form 

sf(s,t,u) =J^ F [^ s (s,t,u) + ^(s,t,u)]+0(p 8 ). (2.1) 

jYf is an overall normalization, the polynomial part J^(s,t,u) contains free parameters describing 
the energy dependence of the processes in analogy to the Dalitz parameters and all the non-analytic 
part of the amplitude connected with the final state interactions is contained in ^(s,t,u). It depends 
on the parameters from the polynomial part and on parameters describing the interaction amplitudes 
for the intermediate states. In the decay region and up to (and including) two loops, the only such 
interactions one has to take into account are the 7T7T rescatterings - the intermediate states other than 
the ones with two pseudo-Goldstone bosons are suppressed to the 0(p s ) order up to a polynomial 
contribution included already in £P(s,t,u). Since we are interested in the amplitude in the decay 
region of the Dalitz plot (i.e. far from the Pn threshold), the contributions of other two meson states, 
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Figure 1: Scheme of the reconstruction procedure. 

like P 71, can be reasonably expanded up to third order in the Mandelstam variables and included 
into &(s,t,u). 

Our reconstruction procedure uses the methods of reconstruction theorem from ^ and will be 
discussed into detail in ph. A brief description of it can be found also in our proceedings [jlj]. From 
the theorem we obtain that the polynomial £P(s, f , u) is a third order polynomial in the Mandelstam 
variables having the same s,t,u symmetries as the amplitude &tf{s,t,u), while the unitarity part 
% (s,t,u) is determined from the single- variable dispersive integrals over the imaginary parts of S 
and P partial waves of all the crossed amplitudes (obtained from the unitarity relation). 

A simple diagram of this procedure is given in Fig. 1. At every order a given parameter 
within the unitarity part (computed from the imaginary part) corresponds exactly to the polynomial 
parameter of the lower order whereas the corresponding parameter within the polynomial part 
contains also contributions (corrections) of the higher order. We can deal with the appearance of 
the parameters of these two orders within the results in two ways: 

• fit order-by-order - we can exactly respect the chiral orders and distinguish contributions of 
every order. We can further simplify the fit by starting with fitting the lower order formula 
and afterwards fit only contributions of the higher orders to the polynomial part. 

• resummed fit - we can include also part of the neglected orders by replacing all the param- 
eters by their highest order values. This possibility is natural in the cases where we want 
to give or keep the physical meaning of the parameters (this will be of use e.g. when fitting 
scattering lengths from the cusp data). 

The amplitudes of all the processes ( \l.l\ > depend on the same %% scattering parameters. We 
can choose them to be the scattering lengths, effective range parameters etc. and write the %% 
scattering amplitude in such parametrization where they keep their physical interpretation up to 
two loops. 
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Figure 2: A comparison between the one-loop chiral perturbation theory result |g] and our one-loop 
parametrization of the amplitude rj — > 3n° on the cut t = u. The real part of our parametrization is fit- 
ted on the chiral results, whereas the imaginary part of it comes from the LO values of our constants. With 
the dashed lines the physical decay region is indicated. 



3. Results 

At present we have a closed analytic form of the parametrization for the Kl^3% and rj —>3n 
decays valid up to two loops (including the pion mass differences). Nevertheless, we present here 
the fit of our results for rj — > 371 decay, calculated in the first order in isospin breaking, where the 
interesting physics appears already there. 

We have put a part of the parametrization of this process already in our previous proceedings 
[j|] and its explicit form will be given in [0]. Let us remind here that at this isospin-breaking order, 
the charged rj — » and the neutral rj — ► 371° decays are related and up to two loops we 

have 6 real parameters describing the energy dependence of this decay together with 4 parameters 
describing KK scattering. We have taken the kk parameters from Stern et al. [§] and fitted the 
results of Bijnens and Ghorbani ||3|] as follows. 

3.1 Order-by-order fit 

Our parametrization with parameters obtained by order-by-order fit should reproduce the chiral 
perturbation theory results ||3|] on the physical decay region. This was indeed verified. On Fig. 2 
and Fig. 3 we plot real and imaginary parts of the 0{p A ) and the 0(p 6 ) amplitudes of rj —> 3k q 
on one particular cut in kinematical variables it = u). The pictures are similar also for other cuts. 
At every order we have taken the constants within the unitary part from the lower order and then 
fitted just the polynomial part. It means that the imaginary part at each order is not fitted to the 
%YT points. The 0(p A ) imaginary part under the nrj (and KK) threshold is reproduced exactly in 
the case when our amplitude is defined consistently to the chiral one (e.g. GMO relations). The 
influence of the nrj and the KK intermediate states is visible already at 0{p A ) above the physical 
region. 

The small shift in the imaginary part of 0(p 6 ) is caused by absorptive contributions from 
sunset diagrams. Their effect can be included by adding a small imaginary part to the 0(p ( ') 
parameters. However, as demonstrated on Fig. 4, to have a good correspondence, we need in fact 
just a small constant shift which is suppressed by two orders of magnitude with respect to the real 
part of the 0(p 6 ) contribution to the constant term. 
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Figure 3: A comparison between the two-loop ^PT result |Qj and our parametrization of the amplitude 
rj — > 3n° on the cut t = u. The real part of our parametrization is fitted on the chiral results, whereas the 
imaginary part of it comes from the values of our constants determined from the 0(p A ) fit from Fig. 2. 
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Figure 4: The same comparison of imaginary parts as in Fig. 3 with adding also a small imaginary part to 
the constant parameter in our parametrization that is fitted from the chiral result. 



3.2 Order-by-order fit vs. resummed fit 

The further illustration can shed light on the sign discrepancy between ^PT ^ and the dis- 
persive determination |Q] of the Dalitz parameter a of the neutral v\ — > 37T° decay The Dalitz 
parametrization of this process is defined according to 

P(. A .)P-aS(i + 2« + ...), z = <3.D 

where Ej is the energy of the i— th pion in the final state. The #PT analysis [||] gives a positive value 
for a whereas the dispersive approaches from [§] (and also experiments) predict negative values. 

Using our order-by-order fit of #PT data (respecting the chiral orders as #PT does), we plot the 
dependence on s and t of the Dalitz plot of \M(s,t, u) | 2 coming from NNLO #PT [|] on Fig. 5 (left). 
To reproduce the dispersive methods [Q], we fit our resummed-fit parameterized \M(s,t,u) | 2 of r] — > 
7r°7r + 7r~ to the chiral 0(p 6 ) results at the Adler zero (s = |M 2 ), where these approaches believe 
the chiral expansion. Using these resumed-fit parameters, we can determine Dalitz parametrization 
for a neutral Dalitz plot on Fig. 5 (right). (This fit should be equivalent to the second numerical 
iteration in the Colangelo et al. [Q] approach - modulo possible different inputs.) As we see on 
Fig. 5 the slope of Dalitz plot (corresponding to a) has really changed the sign. This can also affect 
the computation of the decay widths of these amplitudes and change the value of R. 
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Figure 5: A comparison between our order-by-order fit (left) of tj — > 3ft (corresponding to the ^PT ap- 
proach ^ respecting the chiral symmetry) and the resummed-fit parametrization of this process with the 
parameters obtained by fitting the 0(p 6 ) [[|] rj — * n°n + n~ amplitude at Adler zero (equivalent to the second 
numerical iteration of [||]). 

4. Summary 

We have a general method allowing us to construct a model-independent two-loop parametri- 
zation of the amplitudes of interesting decay modes K — > 3n and v\ — > 371, based only on analyticity, 
unitarity, crossing symmetry, relativistic invariance and chiral power-counting. Our method can be 
implemented in a straightforward way in the presence of isospin breaking (i.e. to describe the cusp). 

We have verified that this description is fully compatible with the previous existing chiral 
two-loop calculation for T] — > 3n in the first order in isospin limit. It is however more general 
and can help to understand also the discrepancy in the sign of slope in Dalitz parametrization 
between the standard chiral calculation and the previous dispersive approaches. It indicates that 
the chiral expansion of the slope parameter a might have a convergence problem and that different 
resummations of higher orders can have a bigger effect on it. This together with the prospects of the 
increase in statistics for T7 — ^ 3?r decays in future experiments ask for a more realistic description 
of the Dalitz plot structure which would go beyond the simple linear + quadratic parametrization 
used so far in the literature. 

References 

[1] M. Zdrahal, K. Kampf, M. Knecht and J. Novotny, PoS (EFT09) 063 [arXiv:0905.4868 [hep-ph]]; 
K. Kampf, M. Knecht, J. Novotny and M. Zdrahal, Nucl. Phys. Proc. Suppl. 186 (2009) 334. 

[2] Crystal Ball/MAMI, Crystal Barrel, KLOE, KTeV, NA48, SND, WASA/COSY, . . . 

[3] J. Bijnens and K. Ghorbani, JHEP 0711 (2007) 030. [See also J. Bijnens, PoS(CD09) 031.] 

[4] J. Kambor, C. Wiesendanger and D. Wyler, Nucl. Phys. B 465 (1996) 215; A. V. Anisovich and 
H. Leutwyler, Phys. Lett. B 375 (1996) 335; G. Colangelo, S. Lanz, E. Passemar, PoS(CD09) 047. 

[5] J. Gasser and H. Leutwyler, Nucl. Phys. B 250 (1985) 539; B. Borasoy and R. Nissler, Eur. Phys. J. A 
26 (2005) 383; M. Bissegger et ai, Phys. Lett. B 659 (2008) 576. 

[6] J. Stern, H. Sazdjian, N. H. Fuchs, PRD 47 (1993) 3814; M. Knecht, B. Moussallam, J. Stern and 
N. H. Fuchs, Nucl. Phys. B 457 (1995) 513; M. Zdrahal and J. Novotny, PRD 78, 1 16016 (2008). 

[7] K. Kampf, M. Knecht, J. Novotny and M. Zdrahal: in preparation. 



6 



